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ABSTRACT 

We report on VLBI, as well as VLA radio observations of the Type Ib/c supernova 2009bb. The 
high radio luminosity of this supernova seems to require relativistic outflow, implying that the early 
radio emission was "engine-driven", that is driven by collimated outflow from a compact object, even 
though no gamma-ray emission was seen. The radio light curve shows a general decline, with a "bump" 
near t = 52 d, seen most prominently at 5 GHz. The lightcurve bump could be either engine- driven, 
or it might represent the turn-on of the normal radio emission from a supernova, driven by interaction 
with the CSM rather than by the engine. We undertook VLBI observations to resolve SN 2009bb's 
relativistic outflow. Our observations constrain the angular outer radius at an age of 85 d to be 
< 0.64 mas, corresponding to < 4xl0 17 cm and an average apparent expansion speed of < 1.74c. 
This result is consistent with the moderately relativistic ejecta speeds implied by the radio luminosity 
and spectrum. 

Subject headings: supernovae: individual (SN2009bb) — radio continuum: general — gamma rays: 
bursts 



1. INTRODUCTION 

Supernova SN 2009bb was discovered by the 
Chilean Automatic S upernova Search Program (CHASE; 
iPignata et all l2009bl laf) on 2009 March 29.9 UT, in the 
nearby spiral galaxy NGC 3278. The radial ve locity of 
NGC 3278 is 2964 km s" 1 (|Paturel et all 12001 . and in 
what follows we will use a round distance of 40 Mpc for 
the galaxy and the supernova. SN 2009bb is in a re- 
gion of high star-formation, approximate ly 4.2 kpc from 
the c e nter of the galaxy in p rojection (jLevesque et al.1 
120101) . iStritzinger et al.l (|2009f) obtained an optical spec- 
trum which showed no evidence for hydrogen, and thus 
SN 2009bb was classified as type Ib/c. The shock break- 
out date is well constr ained to be March 19 ± 1 UT 
(jSoderberg et al.ll2lT10bh . 

Radio emission was discovered using the NRAC0 Very 
Large Array (VLA) on April 5.2 UT, at tim(C3 after 
shock breakout, t of 17 d. The 8.5-GHz flu x den- 
sity was 24.5 ± 1.2 mJy (jSoderberg et al.l 120 10b| ). This 
flux density corresponds to a spectral luminosity of 
^5xl0 28 erg s^ 1 Hz -1 , which is larger than that ob- 
served for any other SN I b/c at a similar time af- 
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ter shock breakout (see ISoderberg et~aIll2010U I2006bl : 
iBerger et al.l [2003). Subsequent VLA observations con- 
firmed the initially high flux density and showed a power- 
law d ecay, with the flux dens ity at 8.5 GHz, Ss.sghz oc 
t~ 1A (So derberg et al.|[2010bl ). Similar decay rates are 
seen for other Type Ib/c SNe, but also for the near- 
est gamma-ray burst, GRB 980425. These radio obser- 
vations gave a position for SN 2009bb of 10 h 31 m 33?87 
and —39° 57' 30 / .'l with an uncertainty of 0.7" in each 
coordinate. We give this position, and all others in this 
paper using J2000 coordinates. 

Radio emission in a supernova is generated by the 
shocks formed as the ejecta interact with the circumstcl- 
lar material (CSM). Radio emission therefore traces the 
fastest ejecta, unlike the optical emission, which traces 
the massive but more slowly-moving bulk of the ejecta. 
Strong radio emission consequently is a sign of particu- 
larly strong interaction with the CSM, which can be due 
either to a particularly dense CSM, such as is seen for 
the radio-luminous Type II SNe, or to particularly strong 
shocks which are caused by relativistic ejecta, whether 
those latter are collimated or not. 

Type I b/c SNe like SN 2009bb are of special inter- 
est because GR Bs have been shown to be associated 
with them (e.g.. iGalama et al.l 119981: iStanek et all 1 2003 1: 
[ Malesani et al1l2004l: iPian et al.ll2006l: iCobb etldlboiol; 
iStarling et al.l 20*To|) . While the optical luminosities of 
Type Ib/c SNe and those associated with GRBs over- 
lap, the GRBs are distinguished by having powerful non- 
thermal "afterglow" emission. In the radio, the after- 
glow typically peaks a few days after the explosion, and 
GRB radio luminosities are observed to be up to a mil- 
lion times higher than thos e of ordinary Type Ib/c 
SNe (jSoderberg et al.l l2006aft . This bright emission is 
the observational manifestation of the substantial en- 
ergy coupled to relativistic ve locities in GRBs . How - 
ever, GRBs are rare events, and lSoderberg et all ([2006b) 
showed that less than 3% of all Type Ib/c SNe have sim- 
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ilarly relativistic outflows. The presence of relativistic 
ejecta, therefore, makes a supernova of particular inter- 
est. The physical mechanism that distinguishes ordinary 
Type I b/c SNe from GRB-SNe remains unknown and 
detailed studies of relativistic Type I b/c SNe are there- 
fore required to make progress. 

In particular, iSoderberg et al.l (|2010b|) showed that 
SN 2009bb's high radio luminosity requires a substan- 
tial relativistic outflow powered by a "central engine" , in 
other words a black hole or a neutron star surrounded 
by an accretion disk which produces collimated outflow. 
The radio spectrum as measured at the VLA and the 
Giant Meterwave Radio Telescope (GMRT) is well fit 
by a synchrotron self-absorption (SSA) spectrum. The 
high luminosity and relatively low turnover frequency 
of ~6 GHz at t = 20 d then imply a b lastwave radius 
of 4.4xl0 16 cm (jSoderberg et alJl2010bt ) and therefore a 
mean apparent expansion speed of 0.85 ± 0.02c, assum- 
ing equipartition of energy between electrons and mag- 
netic fields. Note that these are minimum values for the 
size and expansion velocity, since both deviations from 
equipartition, and the presense of free-free absorption 
(FFA) in addition to SSA would result in larger values. 
No gamma-ray counterpart was detected, but an off-axis 
viewing angle or a low fluence burst cannot be excluded. 
SN 2009bb differs from GRB s in that it occurred in a 
high-metallicity environment (|Levesaue et allfeOlOl ). 

For relatively nearby SNe, a direct measurement of 
the size of the Shockwave is possible with very-long- 
baseline interferometry (VLBI) observations. Such a 
measurement provides a mo del- independent way of mea- 
suring the expansion speed and therefore determining 
the presence or absence of relativistic ejecta, and pos- 
sibly also determining the emission geometry and test- 
ing the assumption of equipartition. Unfortunately, SNe 
which arc both sufficiently nearby and radio-bright to 
allow VLBI imaging are rare events. Relativistic ex- 
pansion was clea rly detected using VLBI in the case 
of GRB 030329 (jTavlor et all [2004). In the case of 
SNe not associated with GRBs, however, the VLBI 
observations so far have confirmed the rarity of rela- 
tivistic ejecta: VLBI observations of two Type I b/c 
SNe which were suspected of having relativistic ejecta, 
SN 2008D and SN 2001em showed only subluminal 
expansion (IBietenholz et all l2009t IParagi et all 2008 ; 



Schinzel et all l2009t IBietenholz fc Bartell I2007L 12005 
Paragi et all 120051 ). In the cas e of SN 2007 g r, rela tivis- 
tic expansion was claim ed bv IParagi et al.l (|2010f >. but 
ISoderberg et al . ( 2010a) showed that a more conservative 
interpretation of a normal, non-relativistic supernova can 
also be reconciled with the VLBI measurements, and pro- 
vides a more natural explanation of the relatively low ra- 
dio and X-ray luminosity. Optical and infra-red spectra 
also suggest a modes t ejected mass and explosion energy 
(jMazzali et al.|[2010^ . whereas relativistic ejecta are usu- 
ally accompanied by large ejected masses and explosion 
energies. 

2. OBSERVATIONS 

We obtained both VLA total-flux-density and VLBI 
imaging observations of SN 2009bb. The VLA observa- 
tions were obtained in two ways. Firstly, by using the 
VLA as a standalone interferometer in parallel with the 
use of the phased VLA as part of our VLBI array on 2009 



June 10 - 11, described below. These observations were 
at 8.4 GHz in the CnB array configuration. Secondly, 
as part of a regular VLA monitoring program for Type 
I b/c supernovae (AS983; PI Soderberg) on 2009 Octo- 
ber 23, at both 8.4 and 5.0 GHz, and with the array in 
the D configuration. The observations were reduced in a 
standard manner, with the flux density scale being set by 
observations of 3C 286, and VCS4 J1036-3744 (hereafter 
J1036-3744; also known as QSO B1034-374) being used 
for phase-referencing. 

The VLBI observations were carried out at 8.4 GHz, 
and lasted for 10 hours with a midpoint of 2009 June 
12.1 UT, or t = 85 d. Our VLBI array consisted of the 
NRAO VLBA (10 x 25-m diameter) and phased VLA 
(130-m equivalent diameter) telescopes, and the Hobart 
(25-m diameter) and Tidbinbilla DSS45 (34-m diameter) 
telescopes of the Australian Long Baseline Array. Un- 
fortunately, due to a bearing failure, the 26-m antenna 
at Hartebeesthoek, South Africa was not available for 
these observations. We recorded a bandwidth of 64 MHz 
in both senses of circular polarization with two-bit sam- 
pling, for a total bit rate 512 Mbit s -1 , with the exception 
of the Tidbinbilla antenna, at which we only recorded 
left circular polarization (IEEE convention). The VLBI 
data were correlated with NRAO's VLBA processor, and 
the analysis carried out with NRAO's Astronomical Im- 
age Processing System (AIPS). The initial flux density 
calibration was done through measurements of the sys- 
tem temperature at each telescope, and then improved 
through sclfcalibration of the reference source. A correc- 
tion was made for the dispersive delay due to the iono- 
sphere using the AIPS task TECOR, although the effect 
at our frequency is not large. 

We phase-referenced our VLBI observations to J1036- 
3744, for which we use the position from the 
Fourth VLBA cal ibrator survey of 10 h 36 m 53M3961, 
-37° 44' 15'.'0662 (jPetrov et al.ll2006h . We used a cy- 
cle time of ~4 min, with ~2.7 min spent on SN 2009bb. 
Due to SN 2009bb's southern declination, it was at rel- 
atively low elevation at the North American antennas. 
The source never exceeded 8° in elevation at the BR, 
HN and NL antennas of the VLBA, therefore we did not 
use the data from these antennas in the final analysis. 
The majority of our remaining visibility measurements 
were made at relatively low elevations, with less than 
25% of our individual visibility measurements being ob- 
tained with both antennas observing at elevations > 15°. 
We show the final u-u-coverage obtained in Figure [U 

For both imaging and model-fitting of the VLBI data 
for SN 2009bb, we reduced the weights of the VLA by a 
factor of 4. As the VLA is more than an order of magni- 
tude more sensitive than any of the remainder of the an- 
tennas, the fraction of visibilities involving the VLA will 
have much higher weight and thus dominate. Although 
reducing the VLA weights incurs a penalty in statistical 
efficiency, it improves the stability of both imaging and 
modelfitting. 

3. RESULTS 

3.1. VLA Total Flux Density 

We describe first the results from the reduction of the 
8.4 GHz VLA interferometric data. Our flux-density un- 
certainties include both statistical standard errors and an 
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Fig. 1. — The n-t)-coverage obtained for our VLBI observating 
run for SN 2009bb on 2009 June 10 - 11 at 8.4 GHz after editing 
and excluding measurements with antenna elevations below 12°. 
The inset shows a detail of the ±10 MA X ±8 MA central region of 
the u-v plane showing the coverage of the shortest baselines. 
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RIGHT ASCENSION (J2000) 

Fig. 2.— An 8.4-GHz VLA image of NGC 3278 on 2009 June 
12. SN 2009bb is indicated. The contours are at -0.1, 0.1, 
0.14, 0.20, 0.28, 0.5, 1.0, and 2.5 mjy bm _1 , and the greyscale 
is labeled in mjy bin -1 . The peak brightness in this sub-image 
was 2640 (j,Jy bm -1 , the background rms 31 fijy bm -1 , and 
the FWHM of the convolving beam, indicated at lower left, was 
3.42" X 2.03" at p a. 14°. The position of the galaxy center 
UPature l ct al. 2003) is indicated by a cross. 

assumed 5% uncertainty in the flux density scale, added 
in quadrature. For our calibrator source, J1036-3744, we 
obtained a flux density of 0.77 ± 0.04 Jy. 

We show the VLA image of NGC 3278 and SN 2009bb 
in Figure [5] The presence of a small amount of extended 
emission due to the galaxy NGC 3278 is apparent. 

By fitting an elliptical Gaussian of the same dimensions 
as the restoring beam and a zero level, we determine the 
8.4-GHz flux density of SN 2009bb at t = 85 d to be 
2.47 ±0.19 mjy, where our uncertainty includes a contri- 



bution due to the uncertainty in estimating the zero-level 
(as well as the aforementioned 5% uncertainty in the flux- 
density scale, all added in quadrature). We estimate the 
contribution from NGC 3278 to be 240 ± 100 fiJy bnr 1 
at our resolution of 3.42" x 2.03" (FWHM). 

For our subsequent VLA observations on 2009 Oct. 23, 
the array was in the lowest-resolution D array configura- 
tion with a synthesized beamwidth of 31" x 6" at 8.4 GHz 
and 47" x 10" at 5.0 GHz. The subtraction of the galaxy 
component of the radio emission was not straightforward. 
We accomplished it by using the image of SN 2009bb and 
NGC 3278 made on 12 June 2009 (Figure©, which had 
adequate resolution and good u-v coverage, as a tem- 
plate. We assume that NGC 3278 does not change with 
time, and that any change in the image between June and 
October is therefore due to SN 2009bb. We convolved the 
template image to the resolution of the new ones, and in 
the case of 5 GHz, scaled the brightness distribution by 
an assumed spectral index of —0.6, and then determined 
the change in the flux density of SN 2009bb. We find the 
flux density of SN 2009bb on 2009 Oct. 23 (t = 218 d) 
was 0.97 ± 0.24 mjy at 8.4 GHz and 1.9 ± 1.0 mjy at 
5.0 GHz. Our uncertainties again include a contribution 
from the galaxy subtraction, as well as those from noise 
and the flux-density calibration, all added in quadrature. 

In Figure[31 we plot the radio light curves of SN 2009bb 
at 8.4 and 5 GHz including the flux density measure- 
m ents described above, as well as earlier values reported 
in lSoderberg et all (|2010b[ ). Note that the VLA flux den- 
sity obtained for the date and frequency of the VLBI ob- 
servations (t = 85 d, 8.4 GHz) shows no noticeable dis- 
crepancy with the remainder of the observed light curve. 

A bump in the light curves is observed near t ~ 52 d. 
It is most prominent in the 5-GHz light curve, which 
shows an increase by almost a factor of 2, reaching a local 
maximum of 12 mjy at t ~ 52 d, which corresponds to a 
spectral luminosity of 2.3 xlO 28 erg s^ 1 Hz -1 . A similar 
although smaller relative increase is seen at 8.4 GHz, 
with the peak occurring perhaps slightly earlier. 

A weighted least-squares fit to all the data up to t — 
220 d shows that on average, the flux density decays 
as follows: S 8 .5GHz oc t- 1A±0A and S'sghz oc t-°- 7±0 - 2 . 
There is a possible flattening of the flux density decay 
at the very latest times (t > 100 d), although due to the 
difficulty of galaxy subtraction the reality of any late- 
time flattening is difficult to confirm. While the overall 
average decay rate is significantly flatter at 5 GHz than 
at 8.4 GHz, this is largely due to the bump being more 
prominent at 5 GHz, with the decay rates being similar 
at both frequencies after the bump. 

3.2. VLBI 

We now discuss the results of the VLBI observations, 
turning first to our phase-calibrator source, J1036-3744, 
as the results obtained for it will inform our interpre- 
tation of the SN 2009bb results. From a deconvolved 
(CLEAN) image of J1036-3744, we recovered a total flux 
density of 687 mjy, with an background rms brightness 
of 1.6 mjy bm -1 . The image is dynamic-range limited, 
rather than limited by thermal noise. The total flux den- 
sity recovered is 10% less than the flux density measured 
at the VLA. It is possible this discrepancy is due to 
10% of the flux density being at spatial scales too large 
for VLBI but too small for the VLA (0.03" - 0.3"), 
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Fig. 3.— The 4.8 and 8.4-GHz radio light curve of SN 2009bb 
from VLA measurements. The measu rements are described in this 
paper or in Sodcrbcre ct al. (2010b). Our uncertainties include 
an estimated 5% uncertainty in the VLA flux-density scale and a 
contribution from the galaxy subtraction. The VLA flux density 
measurement corresponding to our VLBI observations is circled. 
For better visibility, we have shifted the 5-GHz points very slightly 
in time so as to avoid overlapping error bars. We also indicate the 
range of flux densities recovered from the VLBI observations to 
illustrate the discrepancy between them and the VLA flux densities 
(see appendix for details). 

although it is not uncommon to see such discrepancies 
between the flux density scales for VLBI, determined 
from the system-temperature measurements, from the 
more accurate one for VLA measurements, determined 
by observations of calibrator sources such as 3C 286. The 
residual delays and delay-rates we found for J1036-3744 
were moderate, being mostly < 40 nsec and < 5 milli-Hz 
respectively, suggesting that the array was performing 
well, and therefore no particularly strong decorrelation 
is expected. 

We turn now to SN 2009bb. An image made using 
complex weighting with the robustness factor set to —2, 
which is close to uniform weighting, and a u-v taper of 
30% at 120 MA had a total clean flux density of 740 //Jy, a 
peak brightness of 610 /iJy bm _1 , and a background rms 
brightness of 130 fiJy bm _1 . We tried different weighting 
schemes and u-v tapers, but in no case was the total clean 
flux density greater than 800 /iJy. 

For marginally resolved sources, such as SN 2009bb, 
the best values for the source size and VLBI flux den- 
sity come from fitting models directly to the visibility 
data, rather than imaging. We choose as a model the 
projection of an optically-thin spherical shell of uniform 
volume emissivity, with an outer radius of 1.25 x the in- 
ner onj^l. Such a model has been foun d to be appropri- 
ate f or other radio supernovae ( see e.g.. lBietenholz et al.l 
12003 iBartel fc Bietenhoh]|2008j ). For a partially resolved 
source such as SN 2009bb, the exact model geometry is 

11 Our results do not depend significantly on the assumed shell 
thickness, as the effect of reasonable thicknesses different than the 
assumed one is considerably less than our stated uncertainties. 



not critical, and our shell model will give a reasonable 
estimate of the size of any circularly symmetric source, 
with a scaling factor of order unit y dependent on th e 
exact morphology (see discussion in Bartcl ct al. 2002). 
The Fourier transform of this shell model is then fit to 
the visibility measurements by least squares. 

For the interested reader, we give the details of the 
modclfitting results in the appendix. Fitting such a 
model to the strictly phase-referenced visibilities for 
SN 2009bb gives our most accurate estimate of its cen- 
ter position, which is 10 h 31 m 33!8762, -39° 57' 30"022, 
with an estimated uncertainty, dominated by systematic 
contributions, of a few mas. 

The total flux density recovered from the VLBI data, 
whether through imaging or modclfitting, was at most 
^50% of that measured at the VLA. This discrepancy 
might suggest that the source is over-resolved by the 
VLBI observations, or, more precisely, that the missing 
flux density is at angular scales too large to be seen in the 
VLBI observations, but smaller than the VLA resolution. 
The VLBI observations had reasonable u-v coverage even 
for baselines as short as 10 MA (see inset in Figured]), 
and are therefore sensitive to structure up to ^20 mas in 
angular size, while the VLA resolution was ~3". If the 
discrepancy between the VLA and VLBI flux densities is 
to be ascribed to a source resolved in VLBI, then its size 
should therefore be in the range of 20 mas to 3". The 
minimum angular size makes it unlikely that such hypo- 
thetical structure could be related to SN 2009bb, since 
for a supernova age of 85 d and the distance 40 Mpc, 
the implied apparent expansion speed is > 25 c, which 
we consider improbable. 

Another possibility is a source of radio emission unre- 
lated to SN 2009bb but coincidentally so close as to be 
within the VLA beamwidth. We consider also this pos- 
sibility to be unlikely. A flux density of at least 1 mJy is 
required to explain the discrepancy, while the brightest of 
NGC 3278's galactic radio emission seen at the VLA was 
only ~400 /^Jybm -1 , and appears well-resolved (Fig- 
ure [2J . It would seem improbable therefore that there 
would be such a bright and compact unrelated source of 
emission located so close to SN 2009bb. 

The third possibility, which we consider most likely, 
is that there is substantial decorrelation in the VLBI 
measurements, due to uncorrected differences in atmo- 
spheric delay between SN 2009bb and the phase-reference 
sourc43- In other words, the phase-calibration for 
SN 2009bb is of poor quality. Any determination of the 
source size must therefore be interpreted with caution. 
For a sufficiently strong source, the phase calibration 
could be improved by selfcalibration. SN 2009bb, how- 
ever, was not bright enough to allow conventional selfcali- 
bration using an image. As an alternative, we introduced 
the antenna phases as free parameters in the u-v plane 
modclfit. Unfortunately, due to the poor u-v coverage 
and the low signal-to-noise ratio, we were only able to 
obtain an upper limit of 0.64 mas to the source radius 
with this procedure. We found, however, that even with 
this effective phase selfcalibration, we were not able to re- 
cover more than about half of the VLA flux density. Our 

12 We note that inaccuracies in the correlator model could also 
cause coherence loss, but are typically considerably smaller than 
the un-modelled contributions of the atmosphere. 
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TABLE 1 

Size Estimates for SN 2009bb 



Age (days) a 


20 


52 


81 


145 


Radius b (xlO 16 cm) 


4.4 


8.0 


22 


26 


Angular Radius (mas) c 


0.073 


0.13 


0.36 


0.43 


Average Apparent Velocity 


0.85 


0.60 


0.71 


0.69 



a Age in days since shock breakout. 
b Linear size computed from radio spectrum by as- 
suming SSA to be the dominant absorption me cha- 
nism. Radio spectra in[Sodcrbcre ct al. (2010bh. 
c Angular radius (for D = 40 Mpc). 



upper limit on SN 2009bb's angular size is not a formal 
statistical limit because of various assumptions made in 
its derivation, but it is intended as a 3er limit, and should 
be reasonably robust. We describe the exact procedure 
by which we arrive at this limit in the appendix for the 
interested reader. 

Since the age of SN 2009bb at the time of our VLBI 
observations was 85 d, using a distance of 40 Mpc as well 
as assuming circular symmetry, this limit on the angular 
radius allows us to obtain a corresponding upper limit 
on the apparent expansion speed of 1.74 c. 

4. DISCUSSION 

We have made VLA and VLBI observations of the type 
I b/c supernova 2009bb. This supernova was of par- 
ticular interest because the high level of radio emission 
showed tha t substantial material w as ejected at relativis- 
tic speeds (jSoderberg et alj|2010rj) . We discuss first the 
evolution of the radio light curve as revealed by our VLA 
measurements of the total flux density 

4.1. Radio Light Curve Evolution 

Our new data reveal a continuing decline of the radio 
light curve of SN 2009bb at late times (Figure EJ). By as- 
suming that the turnover in the radio spectrum is due to 
SSA, we can estimate the size a t different epochs from th e 
four radio spectra published in lSoderberg et a l. (2010b). 
We give these values of the radius, along with the corre- 
sponding values angular radius and the average expan- 
sion velocity as a fraction of c in Table[TJ These are mini- 
mum radii: if the spectral turnover were due to, for exam- 
ple, FFA rather than SSA, the supernova would be larger. 
However, our VLBI measurement gives a 3cr upper limit 
on the radius of 0.64 mas at t = 85 d, which rules out a 
radius much larger than those derived from SSA. In the 
case of SN 2009bb, much larger sizes and expansion ve- 
lo cities are probably ruled o ut also on energetic grounds, 
as ISoderberg et all (|2010bD showed that the above size 
estimates imply an energy of (1.3± 0.1) x 10 49 erg coupled 
to relativstic ejecta, and much larger amounts of energy 
coupled t o the relat i vistic ejecta are improbable. More 
generally, iChevalierl (|1998[ ) shows that SSA is probably 
the dominant absorption mechanism for Type I b/c SNe. 
In summary, we think it unlikely that the sizes or veloc- 
ities in Table [T] are substantially in error. These radius 
values suggests deceleration between t = 20 d and 52 d, 
but possibly a re-acceleration between t = 50 d and 81 d, 
and a relatively constant speed expansion since then. 

A bump is visible in the light curves, which is most 
prominent at 5 GHz, where it peaks at t ~ 52 d, 
and the flux density increases by a factor of ^2 rel- 



ative to the longer-term decay. In terms of spectral 
luminosity, the bump represents a relative increase of 
-1.2xl0 28 erg s" 1 Hz" 1 . 

What is the origin of this bump? As mentioned, 
SN 2009bb was distinguished by having mildly rel- 
ativistic ejecta. One hypothesis is that the bump 
could be driven by the "engine" , in other words rep- 
resent a renewed injection of energy: similar bumps in 
the lightcurve were seen for GRB 980425 /SN 1998bw 
(jKulkarni et alJ 119981 : iLi fc Chevalier] 119991 ) during the 
relativistic phase, and are thought to be engine- driven. 
The fact that the SSA sizes derived above suggest a re- 
acceleration between i=50 d and 82 d would be consistent 
with this interpretation, in other words that there was a 
renewed energy input at the forward shock somewhere 
around t — 50 d. 

Alternatively, the bump may be related to the colli- 
sion between the bulk of the ejecta, which are almost 
undecelerated, with the leading blast wave, which was 
initially relativistic but decelerates as it sweeps up the 
stellar wind. It is not straightforward to predict the out- 
come of this collision because it essentially depends on 
the poorly known distributions of density and velocity 
in the mildly relativistic, outermost layers of SN ejecta. 
We note, however, that this scenario is similar in many 
respects to that of the bump representing renewed en- 
ergy input by the engine: in both cases the bump in the 
light curve is the result of the collision between the as 
yet undecelerated inner shell and the decelerating exter- 
nal shock. The difference may be in a possible deviation 
from a spherical symmetry, which is expected to be larger 
in the scenario of recurrent jets than in the present case. 
To summarize, we consider the bump in the light curve 
likely to be engine-driven, although the collision of the 
main ejecta with the external shock wave is not excluded. 

Interestingly, even the case that bump is engine-driven, 
i.e., due to recurrent jet activity, the collision of the 
main SN ejecta with the decelerating blast wave is un- 
avoidable. Such a collision might therefore be expected 
in supernovae similar to SN 2009bb or SN 1998bw at 
t = 50 ~ 100 d, and may produce detectable effects in 
the radio lightcurves. 

4.2. Source Size Derived from VLBI Observations 

The primary goal of our VLBI observations was to set 
limits on the expansion speed by measuring the angular 
size of SN 2009bb. Unfortunately, the southern decli- 
nation of the source meant that it was at low elevation 
for most of our VLBI array, which limited the quality of 
the data obtained, and consequently the accuracy of the 
determination of the source's angular size. 

We obtained an estimate of the source outer angular 
radius between and 0.64 mas, suggesting average ap- 
parent expansion speeds between and 1.74c, with the 
limits intended to represent a 3c range (see § !3.2l and the 
appendix for details of our size determination and the 
caveats applying thereto) . This range of apparent expan- 
sion speeds implies a radius at t = 85 d of < 4xl0 17 cm 
a nd a bulk Lorent z factor of < 2.0. 

ISoderberg et al.l (|2010bD showed that deceleration had 
likely occurred relatively early in SN 2009bb, with time of 
the non-relativistic transition £nr — Id, occurring well 
before our VLBI observations. The recent simulations 
of aspherical supernova explosions in stripped-envelope 



6 



Bietenholz et al 



stars by iCouch et all (|2010fl and iKifonidis et"all (|2006l ) 
have also shown that even very aspherical energy releases 
have a largely spherical outer envelope by t ~ 1 d. We 
would therefore expect an approximately spherical blast- 
w ave at the time of our ob servations. 

iSoderberg et al.1 (|2010b| ) also determined that at t = 
20 d, the radius of the supernova was ^4.4xl0 16 cm. If 
we assume a st andard Sedov-von Ne umann- Taylor evo- 
lution (see e.g. JGranot fc Loebll2003T ) after t = 20 d, we 
would expect a radius at the time of our VLBI observa- 
tions (t — 85 d) of ~8xl0 16 cm which is well within our 
observational limits on radius. 

VLBI observations of SNe such as SN 2009bb are cru- 
cial to directly measuring the expansion speeds and per- 
haps the geometry of the radio emission, and thus im- 
portant to confirming a jet model. However, such VLBI 
observations are difficult and generally limited by the 
available sensitivity and u-v coverage, especially if the 
supernova is at a southern declination. They should be- 
come easier in the future with the planned inc reases in 
sensitivity of the VLB A (jUlvestad et alJl2010f ). and by 
future availability of South African a nd Australian SKA 
pathfinder i nstruments, MeerKAT ([Booth et all 120091) 
and ASKAP ([Johnston et al.ll2008fl respectively, as VLBI 
stations. 

APPENDIX 

In this appendix, we describe in detail how we arrived 
at the upper limit on the angular size of SN 2009bb from 
the VLBI observations. We first made an image from the 
strictly phase-referenced VLBI data for SN 2009bb, us- 
ing robust weighting and dropping any data taken with 
telescope elevations below 12°. We found a peak bright- 
ness of only 610 uJy bm _1 for a convolving beam of 
3.4 x 0.8 mas at p.a. 0° (FWHM). Even if we restrict 
ourselves to the baseline-lengths of < 10 MA, we found 
a peak brightness of only 800 ± 160 uJy bm _1 , with a 
convolving beam size of 42 x 10 mas. We note however, 
that due to the uneven u-v coverage, the sidelobe levels 
are quite high, reaching peaks > 70%, and that decon- 
volution is therefore likely not very reliable. 

As an alternative to imaging, avoiding the need 
for deconvolution, we fit a spherical shell model di- 
rectly to the complex visibilities by least squares as de- 
scribed in § 13.21 above. When fitting the strictly phase- 
referenced visibility data, the best-fit model had a total 
flux density of 730 ± 80 \x Jy with an outer angular ra- 
dius, o , of 0.47 mas, and a best- fit center position of 
10 h 31 m 33!8762, -39° 57' 30"022. This is the position 
estimate given in § 13.21 above. 

The total flux density of 730 /zJy recovered from 
strictly phase-referenced data is considerably below the 
value of 2.47 ± 0.19 mJy measured at the VLA (§ l3~Tj) . 
As detailed in § 13.21 we think this discrepancy is unlikely 
to be due to either a supernova which is so large as to 
be over-resolved by the VLBI measurements or to an 
unrelated nearby source of emission. The discrepancy is 
therefore most likely due to decorrelation, in other words 
poor phase calibration. Since most of observations were 
of necessity made at relativly low elevations, difficulties 
in calibration are perhaps not unexpected. 

Our phase-calibration for J1036-3744 is of good qual- 
ity, as attested to by our recovery of 90% of its VLA flux 
density. The delay changes relatively smoothly as a func- 



tion of time, so there is no reason to suspect the inter- 
polation in time between the J1036-374 scans to the in- 
tervening SN 2009bb scans. SN 2009bb, however, is 2.5° 
away on the sky, with the difference being predominately 
in declination. Any un-modeled elevation dependence of 
the delay will therefore result in errors in the delay for 
SN 2009bb, and poor phase-referencing. The sitution 
for SN 2009bb is particularly bad since the difference in 
source position is mostly in declination and since the ob- 
servations are mostly at low elevation where the airmass 
is large. Self calibration in phase would in principle allow 
improving the calibration of the supernova visibilities, al- 
though it does have the drawback of introducing biases 
(see e .g., iMassi & Aaronl 119991: iMarti'-Vidal fc Marcaidi 
20081), which can be severe in the case of low signal-to- 
noise. 

As our signal-to-noise ratio is too low for traditional 
selfcalibration using images, we instead introduce the 
phases of the complex antenna gains as free parameters 
in the model-fitting procedure, using a slightly modified 
version of the AIPS task OMFIT. This procedure has 
the advantage of allowing a more quantitative measure 
of the goodness-of-fit than traditional selfcalibration us- 
ing images. Due to the low signal-to-noise, we fitted for 
only a single phase solution common all 8 intermediate- 
frequency channels. We fix the source position at the 
best-fit position obtained above, and let the antenna 
phases vary on a 30-min timescale. The best fit model for 
SN 2009bb obtained in this way had O = 0.22^1 mas 
and a flux density of 1.34 ± 0.07 mJy (statistical uncer- 
tainties). Even in this case, the total flux density in the 
model was only 54% of that measured at the VLA. 

For a partly resolved source, the fitted source size 
is generally also correlated with the antenna amplitude 
gains. We tested for an additional uncertainty due to 
mis-calibration of the antenna amplitude gains by arti- 
ficially varying individual antenna gains by ±25%, and 
then fitting a model to SN 2009bb as above. The result- 
ing rms variation in o was 0.05 mas. Adding this to the 
above uncertainties in quadrature results in a value for 
o of 0.22±g;g| mas. 

As mentioned, selfcalibration can introduce biases. To 
test for this possibility, we calculated simulated visibili- 
ties from models with various values of 9 and random 
noise at a level corresponding to that in our observa- 
tions. We then fit these simulated visibilities using the 
same procedure as above, including the addition of the 
antenna phases as free parameters. These tests suggests 
that the true value of o is in fact ~18% higher than 
that determined from the fitting. For simplicity, we car- 
ried out this test using a disk, rather than a spherical 
shell model, that the relative bias in 9 Q should be very 
similar. In other words, selfcalibration tends to make the 
source appear more compact. 

Correcting for this bias, we thus arrive at a final, unbi- 
ased estimate of o for SN 2009bb of 0.26±g;?? mas, with 
a 3<7 upper limit of 0.59 mas. However as noted above, 
even with phase-selfcalibration, the fitted flux density is 
still notably below that measured by the VLA, which 
suggests the presence of further decorrelation (or some 
other source of error). As the signal-to-noise is already 
lower than is generally considered safe for selfcalibration, 
using a shorter solution interval is not advisable. 
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As a final test, we fixed the model flux density at 
2.0 mJy, a round value slightly below but near that mea- 
sured by the VLA, and again fitted a model of the source 
as well as the antenna gain-phases. We find that the best 
fit to the VLBI visibilities is not much larger than the 
above estimates, having O = 0.34+q;q4 mas, with a sta- 
tistical 3cr upper limit of 0.64 mas. We note that forcing 
such a large flux density on the model results in a signif- 
icant increase in \ 2 over models with lower flux density. 

The VLBI visibilities, therefore, seem to robustly sug- 
gest a total flux density for SN 2009bb of < 50% of 
that observed at the VLA. Given this inconsistency, the 
biases involved in phase-selfcalibration, and the likeli- 
hood of coherence losses not accounted for by our phase- 
selfcalibration, we suggest a probable range for the outer 
radius of SN 2009bb of < 9 Q < 0.64 mas. The failure 
to recover the total flux density suggests the possibility 
that significant decorrelation remains in the VLBI data, 
but such decorrelation is more likely to increase the ap- 
parent size of the source than decrease it, so our upper 
limit on the angular size should be robust. 

A very similar phenomenon was seen in the case of 
SN 2007gr. Also for this SN, the flux density recovered 
from VLBI observations was considerably lower than the 
total flux density measured by a connected-element inter- 



ferometer. Initially, this discrepancy was interpreted as 
sugg esting a large sourc e size and thus relativistic expan- 
sion (jParagi et al.ll2010h . This was somewhat surprising, 
given that SN 2007gr's peak 8.4-GHz spectral luminos- 
ity was relativ ely low, being ^500 time s lower than that 
of SN 2009bb. iSoderberg et all (|2010a[ ) showed that the 
radio lightcurves and the lack of detectable X-ray emis- 
sion were fully consistent with a normal, non-relativistic 
SN, but were in fact hard to reconcile with relativistic 
expansion. They also re-examined the SN 2007gr VLBI 
data and showed that the low VLBI flux density was ob- 
served on both short and long baselines, and if it was 
to be explained by a large, heavily-resolved source, re- 
quired very large apparent expansion velocities of > 2c. 
They conclude that coherence losses which were larger 
than normal but not improbably so provided an expla- 
nation which as plausible as the original one of modestly 
relativistic expansion for SN 2007gr. In the present case 
of SN 2009bb, some loss of coherence is not unexpected, 
given that the southern declination of the source neces- 
sitated observations made mostly at low elevation. 
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